Nuclear response for the Skyrme effective interaction with zero-range tensor terms 
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I. INTRODUCTION 

Microscopic mean-field approaches are the only ones 
that allow for systematic calculations of binding energy 
and one body observables in the region of the nuclear 
chart that ranges from medium to heavy mass atomic 
nuclei from drip line to drip line [l|. These effective 
approaches rely on a limited number of universal pa- 
rameters, usually fitted on experimental data along with 
properties of infinite nuclear matter derived from realistic 
models Q. 

In the Skyrme-Hartree-Fock formulation, the energy 
of the system takes the form of a functional of local one 
body densities derived from the effective Skyrme inter- 
action ansatz. The commonly used Skyrme interaction 
typically depends on ten parameters and is made of con- 
tact central terms including a spin-orbit interaction. The 
latter, which is controlled by one single parameter, is 
mandatory to obtain the known sequence of magic num- 
bers along the valley of stability. Although this is enough 
to reproduce the global features of nuclei.it was stressed 
from the beginning by Skyrme himself [1, H| that such 
a simple interaction would probably not be sufficient for 
a realistic description of nuclear spectroscopy and that 
a tensor interaction might be needed This last part 
of the effective interaction, made of two contact terms, 
was not considered in most of the early parametrizations 
of the Skyrme force possibly because of the difficulty to 
constrain the corresponding coupling constants. 

In spite of the difficulties related to the adjustment of 
the parameters of the tensor terms, over the years, sev- 
eral attemps have been made for including them. The 
tensor terms as proposed by Skyrme were considered on 
top of the Skyrme SIII Q effective interaction Q or with 
a complete refit of the parameters [1, Q . More recently, 
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the tensor effective interaction has regained some atten- 
tion [h], EH, [H, El, EH EH , partly because it was sup- 
posed to be the key for the reproduction of several spe- 
cific spectroscopic features like, for example, the relative 
shift of proton 137/2 and l/tii/? levels in antimony iso- 
topes In a recent article [13], a systematic study 
of the zero range effective tensor interaction combined 
with a standard Skyrme functional has been made. In 
this work, a set of interactions was built by fixing the 
two parameters of the tensor terms to different values 
while the remaining part of the interactions was fitted 
using the same procedure as for the well known SLy in- 
teraction [ll| EH H(| • It was shown that global features 
of spherical nuclei (masses and radii) and single-particle 
energies cannot be used to clearly mark the boundary of 
a successful domain for the tensor parameters. 

In addition of looking for the parameters that lead to 
the "best fits" to the data, one should also worry about 
the values that could lead to unphysical instabilities of 
nuclei. The tensor interaction energy depending on spins 
and gradients of densities of the interacting nucleons, 
one can intuitively understand that the appearance of 
unphysical finite size domains of polarized nuclear mat- 
ter can be favored for some values of the coupling con- 
stants. Such situations are hard to predict, and to avoid, 
during the fit of the parameters since polarized systems 
break time reversal symmetry while the calculations en- 
tering the standard fitting protocol assume that nuclei 
are spherical and time even. 

A similar kind of instabilities of finite systems was 
encountered and examined in an article devoted to the 
study of effective mass splitting [21j . It was shown that 
the linear response formalism applied to the Skyrme en- 
ergy functional can be used to predict the appearance 
of finite-size instabilities in nuclei. However, only the 
central part of the Skyrme interaction was taken into 
account for the building of the linear response. In the 
present article, we derive the full linear response from 
the energy calculated from a Skyrme effective interac- 
tion that contains spin-orbit and tensor terms. More 
specifically, we always consider that the energy is derived 
from an effective interaction contrary to the spirit of the 
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Skyrmc energy density functional method (Skyrmc EDF) 
for which this link is not required. 

A significant effort is nowadays devoted to the intro- 
duction of new terms in the Skyrme EDF, the tensor term 
being one [l(| [HI, [13 1 other popula r choices being new 
density dependent coup ling s [22. l23l. 0, HH, [26[ or higher 
order derivative terms [271 ] . The number of parameters 
becoming larger, it is of particular importance to define a 
clever fitting strategy Acceptable ranges of variation for 
the parameters can be motivated by linking the Skyrmc 
EDF to realistic interactions using many-body perturba- 
tion theory on top of renormalized low momentum inter- 
actions [H, [2^| . It is also obviously mandatory to locate 
the regions in parameters space that lead to instabilities. 
This is particularly crucial for the spin channels, since the 
corresponding instabilities can only develop when the pa- 
rameters are plugged in calculation codes that allow for 
the breaking of time reversal symmetry 

The linear response function is the tool of choice which 
allows us to avoid areas of instabilities. This application 
will be presented in a forthcoming article, the present 
one being mainly devoted to the derivation of the general 
formulae. 

This work is organized as follows: section [TT] summa- 
rizes the components of the Skyrme interaction which 
includes a zero range tensor, section HTT1 recalls the stan- 
dard formalism of the linear response in nuclear mat- 
ter and discusses its generalization to these new tensor 
terms. In section HV1 we present some numerical calcu- 
lations of responses. Finally, we discuss further possible 
developments in the conclusion. 

II. SKYRME INTERACTION WITH TENSOR 
TERMS 

The usual ansatz for the Skyrme effective interac- 
tion [H, [Til leads to an energy density functional which 



can be written as the sum of a kinetic term, the Skyrmc 
potential energy functional that models the effective 
strong interaction in the particle-hole channel, a pairing 
energy functional, the Coulomb energy functional and 
correction terms to approximately remove the contribu- 
tion from the center of mass motion. The functional 
discussed in this article being applied to infinite nuclear 
matter without pairing, we only consider the kinetic and 
Skyrme potential energy terms. 

Throughout this work, we will use an effective Skyrme 
energy functional, as written in cq. (TAT]) , that corre- 
sponds to an antisymmctrized density-dependent two- 
body vertex in the particle-hole channel of the strong 
interaction. It can be decomposed into a central, spin- 
orbit and tensor contributions 

^Skyrmc = «C + W LS + WT ■ (1) 

We will use the standard density-dependent central 
Skyrme force 

«c(R,r) = t (l + x P a ) 5{r) (2) 

+ Ihil+x^)^' 2 6(r)+S(r)k 2 ] 
+ ta (l + x 2 P CT )k / -5(r)k, 

with the usual shorthand notations for r. R, k, k' and 

Pa [3- 

We will also use the most standard form of the spin- 
orbit interaction 

« LS (r) = i Wo (<ri + <r 2 ) • [k' x 8(t) k] (3) 

which is a special case of the one proposed by Bell and 
Skyrme [H, 130 ] . Finally, the tensor part of the interaction 
is the one proposed by Skyrme 0, [B| 



v T (r) = | t e { [3 («7i ■ k') (<r 2 • k') - {ax ■ <r 2 ) k' 2 ] 8(r) + S(r) [3 (<n • k) (<r 2 ■ k) - (<n • <r 2 ) k 2 ] } 



+ t a 



3 (<7i • k') <5(r) (<7 2 • k) - (<7! ■ <r 2 ) k' • <5(r) k 



(4) 



for which the derivation of the contribution to the to- 
tal energy functional is discussed in detail in refs.[3ll l32j 
as well as in ref.[l7| where the impact of such a tensor 
interaction on the properties of spherical nuclei is inves- 
tigated. 



III. LINEAR RESPONSE FORMALISM 

The linear response function in nuclear matter has al- 
ready been widely developed mainly in the framework of 
Random Phase Approximation (RPA) based on the use 
of an effective interaction [33| . We adopt the presentation 
of the work of Margueron et al. [H| which was devoted 
to the study of the contribution from the spin-orbit term 
to the linear response. 

We consider here the case of infinite matter as a nu- 
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clear medium at zero temperature and unpolarized both 
in spin and isospin spaces. At the mean field level this 
system is described as an ensemble of independent nu- 
cleons moving in a self-consistent mean field generated 
from an effective interaction treated in the Hartree-Fock 
(HF) approximation. For a given density, the momentum 
dependent HF mean field, or self-energy, determines the 
single-particle spectrum e(k) and the Fermi level e(fcp). 

To calculate the response of the medium to an external 
field, it is convenient to introduce the Green function, or 
p-h (particle-hole) propagator G^ a ^(q, w, ki). As it is il- 
lustrated on Figure [TJ ki and k 2 are the initial and final 
hole momenta respectively and q is the transferred mo- 
mentum. We denote by a = (S, M; I, Q) the spin and 
isospin particle-hole channels with S = (S = 1) for 
the non spin-flip (spin-flip) channel, / = (/ = 1) the 
isoscalar (isovector) channel, M and Q being the quan- 
tum numbers related with the projection of the operators 
S and / on the quantification axis. The latter is chosen, 
as usual, as the z axis along the direction of q. 



h - k 2 



q + ki kj 

FIG. 1: Direct and exchange parts of the p-h interaction. 
In the HF approximation, the p-h Green function does 




not depend on the spin-isospin channel (a) and reads |35| 
GffF(g,w,ki) = — ■ — j— jy ■ — p— — — . (5) 



w + e(fci) -e(|ki +q[ 



vquj 



To go beyond the HF approximation one takes into 
account long-range correlations by resumming a class of 
p-h diagrams to obtain the well-known Random Phase 
Approximation [35j . The interaction appearing in the 
RPA is the p-h residual interaction whose matrix clement 
including the exchange part can be written as 



^r'Wr, 



= (q + k 1 ,ki 1 I (a)|V|q + k a ,ki 1 I (a')). 

(6) 

In the general case, the residual interaction is obtained 
by taking the second derivative of the total energy with 
respect to the densities built from the Hartree-Fock solu- 
tions. In the absence of density dependent terms, it can 
also be obtained by standard techniques of particle-hole 
configuration (36|. 

The first important step is thus to determine the ma- 
trix elements for the different parts of the p-h interaction 
from cqs. ©, © and ©. 

A. The particle-hole interaction 



1. Central part of the force 

The central component of the p-h interaction can be 
written in the general form 



(a) 



6 M=0,±1 J 



(7) 



where the w[°^ and coefficients are functions of 

the Skyrme parameters (U,Xi) and of the transferred 
momentum q represented on figure [TJ Their detailed 
expressions have been given by Garcia-Recio et al. HH 
and Navarro et al. [37| for the symmetric nuclear mat- 
ter and pure neutron matter respectively while Hernan- 
dez et al. [3^ gave them for an arbitrary neutron-proton 
asymmetry. The case of symmetric nuclear matter stud- 
ied here is recalled in the Appendix [B] The central part 
of the Skyrme interaction being usually density depen- 
dent, it is not trivialy related to the p-h interaction since 
the coefficients contain rearrangement terms. One 

can note at this level, using only the central part of the 
interaction, that there is no coupling between the differ- 
ent spin and isospin channels. 



2. Spin- orbit part of the force 

To calculate the contribution of the spin-orbit term 
(see eq. ©) to the p-h interaction one has to evaluate the 
matrix element of the spin-orbit interaction. Since this 
term is density-independent there is no rearrangement 
contribution and the result is just adding the following 
term to eq. ([7]) (see Margueron et al. [HI): 



(a., at') 



(<7,ki,k 2 ) = -5 U > w{T) \ —qW 



LS,ph vn-u-i; -jj — \- 1 y g 

xh sl S S ' M [hY^Ch) - k 2 Y^{k 2 ) 
+ SsqS s >xM' \kiY$(ki) - k 2 Y$(k 2 ) 



(8) 
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the factor w(I) = 3 for 1 = and w(I) = 1 for / = 
1 in the case of symmetric nuclear matter. It is clear 
from this expression that the main effect of the spin-orbit 
component is to couple the 5 = and 5=1 channels. 



3. Tensor part of the force 

With the tensor force previously defined (see eq. (HI), 
we have to calculate the antisymmctrized particle-hole 
matrix elements (IQSM;ph\v T \I'Q'S'M';ph). Their 
analytical expressions are summarized in Table [TT] where 
we have adopted the following notation: 



The response function x^Hl^) m the infinite medium 
is related to the p-h Green function by: 



Xrpa(Q, u ) = S I ^~%Grpa(9^M) , (11) 



(2tt) 3 



where the spin-isospin degeneracy factor g is 4 for sym- 
metric nuclear matter. The Lindhard function xhf is 
obtained when the free p-h propagator Grf is used 
in eq. (TTT]). 

Following the notation of refs-JH, [34|, we define for 
any function w, ki): 



= 



k 1 Y^(k 1 )-k 2 Y^ ) (k 2 ) 



(9) 



Even if one can note from Table |TT] that channels with 
different spin projection M arc now coupled in a non 
trivial way. these additional matrix elements are still di- 
agonal in isospin space and act only in the vector channel. 
However, since we include both spin-orbit and tensor in- 
teractions in our approach, it is fundamental to note that 
the tensor component will impact both scalar and vector 
channels via the spin-orbit term. 



B. Response function 

With the particle-hole matrix elements we are now 
in position to solve the RPA problem itself, that is the 
Bethe-Salpeter equation satisfied by the RPA correlated 



Green function G R p A {q,ui,\n\ 



(a') 



d 3 fc 



2 T/ (a,a') 



(2tt) 3 P h 



(10) 



(/(ki)> 



d 3 k! 

(2^)3 



f(kx). 



(12) 



The response function can thus be written in each chan- 
nel (a) as: 

A = 9(G&a)- (13) 
Finally, the quantity of interest is the dynamical struc- 
ture function S^ a \q,uj) which is, at zero temperature, 
proportional to the imaginary part of the response func- 
tion at positive energies: 



' ~ l m Xrpa (l» w ) 



(14) 



C. Response function for the spin-orbit case 

As an introduction to our full calculation we recall here 
some results already obtained by Margueron et ai. [34| . 
When the spin-orbit force alone is included, the response 
function can then be written in the form (using h = c = 
1): 



XHF 

(a) 

Xrpa 



= 1-W^Xo + W^ 



— xo ~ 2k F X2 



xl 



X0X4 



m to 

F 



(15) 



xl 



6ir 2 kF 



q 2 xo 



3tt 2 



w. 



(n) 



Xo • 



In this expression k F is the Fermi momentum whereas 
m* denotes the effective mass of the nucleons. The func- 
tions xo: Xi an d X4 are generalized free response func- 
tions, defined in ref.[33| and written in Appendix [Dj sec 
eq. dDTJ, and Xhf = 9Xo- 



The explicit expressions of the W{ coefficients are 
given in Appendix [C] where the coupling between the 
5 = and 5=1 channels induced by the spin-orbit 
interaction can also be clearly seen. It can be noted that 

the coefficients w[ a ^ are now complex functions of q and 
uj since different moments xn enter their expressions. If 
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we replace W{ in cq. (| 16[) by Wy we obtain the results 
of ref . [33| , related to the central part of the interaction, 
as it should be. 

D. Response function with the tensor part 

As already mentioned, the tensor interaction couples 
vector channels with different spin projection while the 
spin-orbit interaction couples the scalar (5 = 0) and vec- 
tor (5 = 1) ones. The consequence is that wc obtain 
a non trivial system of coupled equations for the RPA 
problem. As an illustration let us consider explicitly the 
case of isospin 1 = and (5, M) = (1,1). In that partic- 
ular channel, the Bethe-Salpeter equation for G R p A (ki) 
(we omit the isospin index and dependence on q and us for 
sake of simplicity) exhibits terms which typically read: 

. G^A* 

• G ffF (fc 1 )(fc 2 G< ? 1 P 1 ] 1 ), 

• fc 1 yj 1) *G ffJ ,(fc 1 )( kY^G^'p } A ) , 

. G ffF (fc 1 )(fc 2 r ( 1 1) r ( 1 1) G( ? 1 i ;i ) )- 

Thus, the determination of ( G R p A ) requires the knowl- 
edge of some other unknown quantities. This leads to a 
large system of coupled equations for: 

• \ U RPA I ' 

• ( k 2 G RPA ) , 

• <kY (1) G {1 - M) ) 



. (fc 2 |y i (1) | 2 G^ ) )forM = -l,0,l, 

and 

• ( kY 1 G R ' P X A ) , 

• ( kY 1 G R p A ) , 

• ( kY_ x G RPA ) , 

that is 21 unknown quantities. Fortunately, since the 
multipolc expansion of Ghf only implies terms with 
Yq , the integration over k\ cancels all terms of the form 
(f(k)n {M ,M>,... } Y£ ) Y$...G HF )withM+M'+...^0. 
Moreover some unknown quantities can be expressed 
through the others and we can reduce the 21 coupled 
equations to three systems of four coupled equations for 
the following variables: 

X^RPA It '-'RPA I ) \ KI '-'RPA h 

\ K l r I '-'RPA 1 ~ \ K \ 1 l I U RPA )i 

for each value of M. The calculations are straightforward 
but tedious. The same procedure has to be repeated for 
isospin I = 1 for 5 = 1. Indeed, only the channels with 
5 = are less involved. 

The results are now quoted for each spin-isospin chan- 
nel in a form that exhibits the symmetry properties ap- 
pearing in Table HT1 



• For 5 = 0, 1 = 0: 



XHF 
(0,0) 
Xrpa 



~(o,o) 



2xo 



W, 



(o,o) / q 



-Xo - 2k F X2 + 



(0.0) 



Wr 



(0,0) 



W^ 0) k 2 F 



X 2 - XoX4 + 



kp 



xl 



67T 2 fc, 



-q 2 xo 



(16) 



• For 5 = 0,7= 1: 



XHF 

(0,1) 
Xrpa 



~(o,i) 
W t xo 



2xo 



W, 



(04) / ff 



m*k 3 t 



-Xo - 2k F X2 



(0,1) 



Wr 



(0,1) 



W^k F 



xl 



X0X4 



fe p 



Xo 



67T 2 fc, 



-g 2 xo 



(17) 
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For S = 1,1 = 0,M = ±1: 



XHF 
(1,0,±1) 
Xrpa 



l + (t e -5t ) 



ir. 



(1,0) 



2tt 2 

(te - 5t ) 
W^° ] - (t e - 5t ) 



2 ~(1.0,±1) 

- H 7 ! Xo 



9 2 

yXo 



1 + (t e - 5t ) 



m*kp 



- 2k\ X 2 + (te - 5t )^-^ (XO - X2) 



X2 - X0X4 + 



Xo 



6n 2 k j 



-q Xo 



+ 2 



W 2 (1 ' 0) + 2(t e - 5t„) + X^ ^ 1 ) 



9/1- 



3tt 2 



W 2 (1 ' 0) + 2(t e - 5t„) - X(i,°,±i) 



2m*wX< 1 '°' ±1 > 



2m*k 3 F 
3tt 2 



^(1,0,±1) 



O TllUJ 

2 + ^r 



Xo 



3tt 2 



(18) 



W 2 (1,0) + 2(t e - 5t ) - X(i,°,±i) 



• For S = 1,1= 0,M = 0: 



XHF 

(1,0,0) 

Xrpa 



1 2 (^e ^t ) 



2tt 2 



W^ ] + \{t e -U ) 



~(1,0,0) 

- W x xo 

1 - (t e - 5f„) 



9 2 

yXo 



m*k F 
2ir 2 



2fcf,X2 - (f e - 5t )^-^ (xo - X2) 



X 2 - XoX4 + 



Xo 



67T 2 fcj 



-q xo 



iy 2 (1 ' 0) -(t e -5t Q ) + x( 1 ^°) 

1 - ^[W^' 0) - (t e - 5t G ) - Xd,",")] 



Yo 



Sir 2 



~ X (1,0,0) 


2 


q m*to 


') 






[2 q 





Xo 



(19) 



1 _ !^[M/ 2 (1 ' 0) - (t e - 5t c ) - X(i,o,o)] 



• For S = 1,1= 1,M= ±1: 



XHF 

(1,1, ±1) 
Xrpa 



~ / s TTh k rp 



2 ~(i,i,±i) 
- Xo 



IT". 



(1,1) 



3(te - to) 



-Xo 



l + 3(t e -t ) 



24X2 + 3(t e - to) 7 ^- ( Xo " X2) 



+ 



W^ 1] - 3(t e - to) 



kp 



Xi - XoX4 + 



kp 



Xo 



67T 2 fcj? 



<Z Xo 



+ 2 



iy 2 (1 ' 1) + 6(te-t ) + x( 1 - 1 ' ±1 ) 

1 - n 4^[W^ 1] + 6(t e - t ) - Xd,i,±D] 



Xo 



q J 3-7T 2 



x (i,i,±i) 



1 _ ^[W^ 1] + 6(t e - to) - Xd,i,±D] 



(20) 



7 



• For S = 1,1 = 1,M = 0: 



XHF 
(1,1,0) 

Xrpa 



l-f(f e - to) 



^2 (1,1) + |(ie-<o) 



m*fcr 



2tt 2 



~(i,i,o) 
W 1 xo 



q 2 

yXo 



l-3(t e -*o) 



m*k 3 1 



2tt 2 



m*fc? 



- 2fc^ X 2 - §(i e - *»)— r 1 - (Xo - Xa) 



W 2 (1,1) + |(t e -*o) 



pi 



X 2 - XoX4 + 



Xo 



67T 2 fcj 



-9 Xo 



W$ L ' 1) -3(t tl -to)+Xl 1 < 1 > ') 



1 - 



TO*H r T Tr(l,l) 



- 3(te - to) - XC 1 ' 1 ' )] 



\0 



. m~w \ m*kp 



X(i,i,o) 



Xo 



[W 2 U ' iJ - 3(t e - to) - X(M,0)] 



(21) 



~{SJ) ~(S,I,M) 

The coefBcients for S = or W 1 for 5=1 are defined as 



;(o,o) 



;(o,i) 



:(i,o,±i) 



;(i,o,o) 



;(1,1,±1) 



W 1 



;(l,i,0) 



ft ~ ft 



1 + g 2 ^ - ft)(^ 2 (1 ' 0) - (7t e + 13t„)) 



^(l,0,±l )+iri+ ^ 2 _ t 2 



i + g 2 (ft-ft)[W 2 (M) +3(i e -f )] 

24m* fcl 



5tt 2 



+ 9(i e + 3i )V(A> - 2/3 8 + ft) + 2X^ ^ 



q m u) 
2 + — 



W^fi) _ i Kl + i K2 _ tl^^t + 2X (i,o,o) 
^^) + 3 ^ + 9^-(t e -to) 2 ^ 



+ 9(i e - i )V(ft - 2/3 8 + r ) + 2X( 1 - 1 ' ±1 ) 



q m u> 
2 



(1,1,0) 3 17/ ,91// 



12ra*£;| 
5tt 2 



+ 2X( 1 ' 1 '°) 



(22) 
(23) 

(24) 
(25) 

(26) 
(27) 



The functions Pi are given in Appendix iDl and the coefficients Ki, K[ and X^ > r ' M ' are 



K\ = 2(t e + 5t )q 2 + 6 



K' 



K 2 = -(t e -5t ) 



2 I l l: kp 



9, 9 3 9 3 / to*w 

— ki + -q 

10 F 8 q 2 \ q 



2{t e + t )q 2 + 6 ( — ) (t e -t„), 



(te-U 2 ^ 



(28) 
(29) 
(30) 
(31) 
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and 



x (i,o,±i) 
x (i,o,o) 

x (l,l,±l) 

x (i,i,o) 



2„2 



3Gf a q 



02 -fc 



l + q*{(3 2 -fo)(w£' 0) + \{t e -ht o )) 



2„2 



= 9(t e -t y q 



1 + g 2 (/? 2 - /3 3 )(IF 2 (1 ' 0) + 5t e + 23i D ) 

02-03 



i + g 2 (0 2 - 3 )(iF 2 (M) + §(t e -i o )) 



2„2 



18(te -U g 



'2 ' 2* 



l + g 2 (/? 2 -/3 3 )(IF 2 (M) -9(i e -t )) 
I 



(32) 
(33) 
(34) 
(35) 



It is important to note that the above response func- 
tions have exactly the same structure as in HF, with or 
without the spin-orbit interaction. Moreover we can see 
that the response function depends on different linear 
combinations of the parameters t e and t a that will lead 
to non trivial effects. 



IV. DYNAMICAL STRUCTURE FUNCTION 
WITH TENSOR CONTRIBUTION 

As an example for the effect of a zero-range tensor force 
in the p-h interaction, we have calculated the nuclear 
responses in (S, I, M) channels of symmetric infinite nu- 
clear matter for parametrization T44 of the Skyrme inter- 
action built by Lesinski et al. (l7j . The parameters of this 
force are given in Table HI The dynamical structure func- 
tions S^^q, uj) (defined in eq. (fl4]l ) calculated for q = kF 
and at the saturation density p = p sa- t = 0.16 fm -3 are 
shown on Figure [2j To clearly isolate the effect of the 
tensor part of the force, the functions S^ a ' (q, u>) are plot- 
ted for two cases: a) with no tensor contribution but in- 
cluding the spin-orbit part; b) with the full force. The 
first case allows to compare with the previous results of 
Margueron et al. f34j but are presented here for the four 
channels of the symmetric nuclear matter. 



to 


x 


ti 


Xl 


ta 


X2 


-2485.67 


0.721557 


494.477 


-0.661848 


-337.961 


-0.803184 


ta 




7 


Wo 


t e 


to 


13794.7 


1.175908 


1/6 


161.367 


173.661 


7.17383 



TABLE I: Parameters of the interaction T44 as denned 
in equations ([2]), Q and (UJ). 

In S = channel, the tensor terms do not affect qual- 
itatively the response. Strictly speaking, this comment 
only applies to the case of the T44 tensor parametriza- 
tion but several tests performed using other Tij tensor 
interactions discussed in ref.JlTj exhibit the same quali- 
tative behaviour. The situation is quite different in S = 1 
channels: the effect from the tensor terms is large what- 
ever the value of the spin projection M is. Actually, 
depending on the values of the transferred momentum q 




tfl 0.0015 



> 0.001 
Ej 0.0005 



B— 1, 1— r 


1 

— M=0, (a) " 

— M=l, (a) ■ 
■ ■ ■ ■ M=0, (b) " 

— - m=i, (b) ; 





S= 1,1=0' 




M=0, (a) 
M=l, (a) 
M=0, (b) 
M=l, (b) 



100 

CO [MeV] 

FIG. 2: Dynamical structure functions (q, uj) calculated 
for the Skyrme tensor parametrization T44 [171 ] . The upper 
(middle and lower respectively) part of the Figure concerns 
the channels S = 0, S = 1 and I = 0, S = 1 and 7 = 0. 
For each (S,I,M) channel the responses (in MeV -1 fm~ 3 ) 
are plotted as a function of u> (in MeV). All the responses are 
calculated for q = fcj? at saturation density. Cases (a) and (b) 
are discussed in the text. 



and the density p, the response functions increase signifi- 
cantly and diverge at finite q for a certain critical density 
p c . This divergence reveals the presence of instabilities 
observed in nuclei 21(, with the appearance of domains 
with typical size of the order of 2Tr/q. Even if a one-to-one 
correspondence between infinite matter and finite nuclei 
is obviously not correct, it remains that the center of a 
nucleus explores, because of fluctuations, not only the 
saturation density but also some larger values for which 
one may observe a divergence of the response functions, 
and then, possibly, the appearance of finite-size instabil- 
ities in the nucleus. 
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V. SUMMARY AND CONCLUSIONS 

We have derived the contribution from the tensor 
terms in the Skyrme effective interaction to the RPA re- 
sponse function. We have shown that the formal struc- 
ture of the response function is the same as without ten- 
sor terms, although with the latter, all channels are cou- 
pled in a non trivial way. The simple example presented 
here, using the interaction T44, shows that the effects of 
the tensor contributions are strong in vector channels. 

We have shown that the dynamical structure functions 
S( a '(q 7 uj) becomes large for finite values of q. This indi- 
cates the vicinity of a pole related with a finite size insta- 
bility for given values of q and p in infinite matter and, 
possibly, in finite nuclei. A systematic study of the criti- 
cal densities is in progress in order to determine if the link 
between the divergences of xb,pa and the instabilities en- 
countered in nuclei at the Hartree-Fock approximation is 
robust. 

Another important point under study is the identifi- 
cation, directly from the Skyrme energy functional, of 
the origin of each tensorial contribution in the response 
functions. In the same spirit, a detailed study of sum 
rules can enlighten the contribution of the tensor for var- 
ious physical situations (see (39|). Finally, applications 
to pure neutron matter can be very important (see by in- 

stance H m m m, 0, m H S3, 5a Hi ML HI m ) • 

However, above formulae arc no longer directly usable 
and have to be adapted to that specific case. Work in 
that direction is also in progress. 
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APPENDIX A: COUPLING CONSTANTS OF 
THE SKYRME ENERGY FUNCTIONAL 



The Skyrme energy density functional used in this ar- 
ticle is a functional of local densities and currents (with 
q = n, p for neutrons and protons), i.e. the particle 
densities p q (r), the kinetic densities r 9 (r), the current 
(vector) densities j g (r), the spin (pseudovector) densi- 
ties s g (r), the spin kinetic (pseudovector) densities T g (r), 
the spin-current (pseudotensor) densities J q ^ v (r) and the 
tensor-kinetic (pseudovector) densities F g (r) that have 
been defined in Lesinski et al. [l7l , where the definition of 
the vector spin current density 3^ (r) = J(r) was also re- 
called. We will use also the isoscalar and isovector densi- 
ties defined from the proton and neutron densities respec- 
tively as: p (r) = p n (r) + p p (r) and pi(r) = p„(r)-p p (r), 
and similar for all other densities. 

The Skyrme energy functional representing central, 
tensor and spin-orbit contributions is given by: 



kyrmc 



£c+£: 



LS 



£ T 



/d 3 r^{ C P M Pt + C^pt^Pt + CI ( Pt r t - j 
J t=o,i ^ 



C°[po] s 2 t + C7 S (V • s t ) z + Cf ' s t • As t 



(Al) 



■a 



VJ 



[ptV-Jt + st- (V xj t )]| 



c* 



si 



2 Jt,H"Jt,vp 



All the coupling constants can be expressed as func- 
tion of the parameters of the Skyrme interaction given 
in eqs. ([2]), ([3]) and ([4]). Some of the coupling con- 
stants are fully defined by the standard central part of 



the Skyrme force: Cf = A p t , Cf p 
and Cf 



/ t — Af , or by its spin-orbit part: C t VJ = Af J , 
Two coupling constants only depend on the tensor part 



of the interaction: C t Vs = BY S , and C, 



Bf. Fi- 
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nally, two coupling constants are the sum of contributions metric nuclear matter (p n 
from both central and tensor forces: CT — AT + BT and 



ip and pi = 0): 



C As = A As + B A.^ 

The coupling constants of the Skyrme EDF which come 
from the central and spin-orbit parts of the interaction 
are given in terms of the parameters by: 



o 

VJ 



- |*0 + 



" 64 



•'•() 

1 



"4*0(2 



2_i (l 

24 1,3 V 2 



Al = 

A\ = 

K = 

At = 

aAs 

A Q — 

Af S = 

Al = 



■^3) Po( r ), 

|<l(3 + *l) + 55*2(| + X 2 ) , 
L + i*2 (|+aj 2 ), 
-*i) + |f 2 C 



3 i 

16*1 



it 



'o' t 



32 11 
64*1 + 



\+X X ) + ft 2 (± + X 2 ) 
n(± - X ) - s*3(| ~ ^3) Po( r ) ) 
- i*3 Po( T ) 



64 "-2 : 



' 16*1 + 16 12 

4 Wb. 



3 - X X 



52*2(2+^2), 



4 " 1 



(0,0) 



(1,0) 



1^(0,0) 



4^2 



(0,1) 



1^2 



(1,0) 



W2 



(1,1) 



2<° + A^( 7 + l)( 7 + 2) A 



(0,1) _ 



2 /1 o 



2Af + 2A^p2 
2Af" + ±A[] 9 2 . 

+ ^0 ] ? 2 : 



(1,1) _ 



2Af - 
[2Af s 

^0 ! 
^1 ! 

A T 

^0 ! 



APPENDIX C: W (a) COEFFICIENTS: 
SPIN-ORBIT PART OF THE FORCE. 



Due to the density dependence of some coupling con- 
stants it is also useful to define the following coefficients 
which occur in the definition of the coefficients [l|: 



When the spin-orbit force is taken into account in the 
RPA formalism, one has to define with the help of the $ 
moments [H[ given in the Appendix [Dl 



A? 



A 



8,0 



A p ^pl, 



• for S = 0, the coefficients W{ 



(S,I). 



In each case where we have considered a single density 
dependent term, the generalization to more than one den- 
sity dependent term is straightforward just adding new 
density dependent terms in the corresponding coupling 
constants. 

Finally, the coupling constants of the Skyrme EDF 
which come from the tensor part of the interaction are 
given by (Table I in [12]): 



(0,0) 



(0,1) 



(0,0) 



2„4 



l + wf-VCfc-ft) 



0a - A 



for S = 1, the coefficients W{ 



(S,I,M). 



3* ), 
3t„), 



■ 8 (*e 
g(*e 
32 (* e *o) : 

^-(t -t ) 

32 \ £ ° ) ' 



Bf 
Bl 
B^ 

bY s 



8 (* e 



a(*e 



!2 (^*< 



32 



3,2 



to), 
to), 
+ t a ), 
+ t a ). 



APPENDIX B: W a COEFFICIENTS: CENTRAL 
PART OF THE FORCE. 

With the coupling constants defined in the Skyrme en- 
ergy density functional and recalled in Appendix lAl the 
W^ a ' coefficients take the following expressions for sym- 



(1,0,1) 



(1,0,0) 



(1,0,-1) 



(1,1,1) 



Wi 



(1,1,0) 



(1,1,-1) 



1,0) 



(1,0) 



(1,0) 



(1,1) 



ir; 



(1,1) 



(1,1) 



?wV 

2 °% + WfV(/32-03)' 



V n V 

IWnV 

2 °% + wf'V(/32-/?3)' 



W V 



02 -/?3 



l + WfVCA-A) 
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APPENDIX D: GENERALIZED LINDHARD 
FUNCTIONS 

Following ref. [33[ , the definition of the generalized free 
response functions is: 



The explicit expression can be found in ref.[33j for i 
0. 1, 2. With X2 and X4 we can compute the different 
moments of Ghf occurring in the calculation. As Garcia- 
Recio et al. (33j , we introduce the functions Pi as: 



X2i(q) 



2(2tt) 4 



4 

k],< 



(Q3 



k 2 



Ghf(<73,9) ■ 
(Dl) 



[A, i = 0, 



d 4 q 3 
(2tt)* 



G 



HF 



1 q qs q| (q ■ qs) 2 (q ■ qshl qf (q-qs) 3 (q ■ qs) 4 (q ■ q3) 2 q 2 



r 



with the functions pi written as 



0o 
2k/3 1 
4k 2 2 



4fc 2 /3 3 



8ft 3 /3 4 
16ft 4 /? 5 



16fc 4 /3 7 



16fc 4 /3 8 



Xo , 

(f - fc )xo , 



m*kF 
6tt 2 



2kv(k - v)xa + (y — k)\2 + - „ ft . 



3tt 2 



X4 - 4fc^X2 



3 /3 6 = (1/ - ft) 3 X o + (3ft - 1/) 



m*kF 
6tt 2 



m*kp 



2tt 2 



fc 2 + I + I(2fc-,) 2 



(1/ - fc) 2 X 2 - - fc) 2 Xo 



6tt 2 



[1 + 2ft(3ft - z/)] 



where ft = of- and 1/ = 

For completeness, we now quote the different moments 



of Ghf encountered in the RPA-cquations: 

( Ghf } = 0o > 
(ftF (1) G^) = (Z JJ/3 1 , 



(fc 2 GW> = g 2 /3 2 



> 2 
'ft 2 



y, 



(i) 



(i; 



G H F)=q 2 --P 3 , 
An 

' G HF ) = q 2 ^-{P2~fo), 



(fc 3 y (1) G^> = ^/^:/34, 

3/2 



[k 3 
[k" 



Y 



(i) 



y (1) Ghf ) = q 3 [i L ) As , 



4tt 



(i) 



Y (1) G HF ) = q 3 \l '-Z-f (Pi- 06) 



3 3 



47T 87T 



(ft 4 G HF ) = q i 0-, 



Y 



(i; 



(ft 4 



(1) 



(1) 



(ft 4 



(1) 



G//F) = g 4 ^(/?8-/3 7 ), 

G H F) = q 4 -^(0 5 -2P 8 + P 7 ). 
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TABLE II: Particle-hole matrix elements for the tensor part of the effective force. Only the vector channels are concerned in this case and a factor 2SsiS S i 1 8 ii i is 
implicit for each term. 
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